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Experimental Verification of Theory of Landing Impact 

By Walter Ramberg and Albert E. McPherson 



Drop tests of an idealized wing and alighting gear were made to provide an experi- 
mental check on methods for computing the transient bending stresses in the wing produced 
by a symmetrical vertical landing impact. The model was dropped in a nearly strain-free 
condition to make contact at a point below the center of gravity. The forces in the alight- 
ing gear, acceleration at the "fuselage", and bending strains in the wing were recorded as a 
function of time. 

According to the statistical theory of Biot and Bisplinghoff , the computed maximum bend- 
ing stresses, using the first three flexural modes, were found to be 43 to 137 percent greater 
than the measured values. Using the actual forcing function reduced the difference to less 
than 20 percent. 



I. Introduction 

The tests described in this paper ] were made 
for the Bureau of Aeronautics, United States 
Navy Department, to provide an experimental 
verification of analytical methods for determining 
the transient oscillations in the structure of an 
airplane during landing impact. 

The practical importance of this problem for 
the safe operation of large airplanes has been 
stressed by Biot and Bisplinghoff [1] 2 , Keller [2J, 
and Yorgiadis [3]. 

The analysis of the transients during landing 
impact is complicated by the fact that these 
transients involve many natural modes of vibra- 
tion of the airplane, and response in each mode 
depends on the force-time curve at the point of 
contact. The force-time curve will vary from 
one landing to the next of a given airplane. In 
view of these complications, Biot and Bispling- 
hoff proposed an ingenious statistical approach to 
the landing problem. In this approach the vibra- 
tion of the structure in a given mode is reduced to 
that of an equivalent linear oscillator, and the 
maximum amplitude in that mode is estimated 
from an envelope of " dynamic response factors", 
which bounds the response to impact force-time 

1 This paper was presented before the Sixth International Congress for 
Applied Mechanics in Paris on Sept. 27, 1946. 

2 Figures in brackets indicate the literature references at the end of this 
paper. 



curves of any shape that may be expected in the 
landing. An upper limit to the resultant ampli- 
tude is obtained by adding up the maximum 
amplitudes in the various modes. 

Application of Biot and Bisplinglioff's statistical 
approach involves the following assumptions, 
which may affect the accuracy of the result: 

1. The maximum amplitude in the various 
modes is added up without regard to phase differ- 
ences. This will lead to a resultant that may be 
considerably larger than the resultant when phase 
differences are taken into account. 

2. The most severe impact force-time curve 
during the landing approaches in effect one of 
the impact force-time curves that were used by 
Biot and Bisplinghoff to derive their envelope of 
" dynamic response factors". Further measure- 
ments of landing impacts in service and of impact 
force-time curves in drop tests of landing gear 
may lead to modifications of the envelope of 
dynamic response factor. The envelope should 
be raised if service showed more severe impact 
force-time curves than those assumed by Biot 
and Bisplinghoff. It could be lowered if they 
were found to be consistently less severe. 

3. It is sufficient to confine the analysis to the 
first few modes of vibration. Inclusion of modes 
of higher order than about the fifth is impracti- 
cable, because they are either unknown or too 
difficult to determine. This may lead to an exces- 
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sively low value of the resultant if the higher 
modes contribute a significant proportion to the 
resultant. 

4. The force-time curve at the landing gear is 
independent of the flexibility of the airplane 
structure. Actually there may be an appreciable 
coupling between the elastic deflection of the 
structure and the action of the landing gear. 

5. The effect of damping is negligible. Damp- 
ing may introduce coupling between the various 
modes of vibration of the airplane, thereby making 
it impossible to analyze the vibration in any given 
mode apart from that in all other modes. 

6. Nonlinear effects are negligible. Nonlinear 
effects caused by buckling of the sheet, slipping 
of rivets, exceeding of the proportional limit of 
the materials in portions of the structure may 
introduce coupling between various modes. 



II. Description of Model 
1. Wing 

The wing was designed to have a mass and 
flexural-rigidity distribution approximating that of 
a Xo-scale model of a large military airplane. The 
assembled model is shown in figure 1. The wing 
was a tapered box beam of rectangular section. It 
was constructed of aluminum alloy sheet and 
angles fastened with %-in. rivets. The nominal 
cross-sectional dimensions at several stations along 
the wing are shown in figure 2. 

The flexural rigidity, EI, at various stations 
along the completed wing was computed from 
measurements of extreme fiber strains when the 
wing was subjected to a known bending moment. 
The results are given in figure 3. 

The mass distribution of the wing was measured 




Figure 1. Idealized model. 



An experimental verification of Biot and Bis- 
plinghoff's analysis was decided upon in order to 
determine the adequacy of the assumptions and 
to indicate the minimum number of modes of 
vibration that must be included to estimate 
maximum bending moments and accelerations at 
various stations along the wing. 

The experimental verification was started with 
the tests described in this paper. These tests are 
concerned with the simplest case, namely that of 
measuring the flexural transients in a symmetri- 
cally tapered model wing when the model is sub- 
jected to a vertical impact force directly below 
the center of gravity. Four engine masses were 
mounted symmetrically on the wing so as to sub- 
ject it to bending without torsion during the 
landing impact. 



by a volumetric method. This involved stripping 
the wing of the engine weights and measuring the 
change in weight as the wing was lowered into a 
large container of water. The measured values of 
mass per inch are shown in figure 4, a. Figure 4, b 
shows the magnitude and location of concentrated 
masses corresponding to the engines and fuselage. 

2. Alighting Gear 
(a) Damper 

Adjustment of the damping force over a wide 
range was made possible with a fluid damper built 
into the center of the model. The damper, 
figures 5 and 6, consisted essentially of a piston 
forcing fluid out of 8 adjustable ports in a cylinder. 
The ports were located in the side of the cylinder 
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Figure 2. Nominal cross-sectional dimensions of model wing. 
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Figure 3. Flexural rigidity distribution for model wing. 
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Figure 5. Schematic diagram cf alighting gear. 

wall at 0.07-in. intervals along the path of travel 
of the piston. As the piston traveled up the cyl- 
inder, it successively closed off the ports. The 
shape of the curve of damping force versus piston 
displacement could be changed by adjusting the 
initial opening of the various ports. 

(b) Spring 

The cushioning action of the air pressure in an 
oleo strut was approximated in the model by a 
pair of double cantilever springs shown near the 
bottom of figure 5. The springs were designed to 
have a wide range of adjustment of stiffness. The 
stiffness was changed by moving the length-ad- 
justing block either closer or farther from the 
vertical centerline of the model. The springs were 
attached between the landing foot and the root of 
the model. The initial load applied by the spring 
between root and wing could be adjusted with the 
spring-loading bolts shown in figure 5. 

(c) Landing Foot 

The landing foot of the model is shown at the 
bottom of figures 5 and 6. It had a conical impact 
surface, which was intended to give a nonlinear 
force-displacement relation as the foot com- 
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pressed the rubber-surfaced landing pedestal. A 
nonlinear force-displacement relation was desired 
to simulate the nonlinear characteristics of the 
tire in an actual airplane. The shape of the force- 
displacement curve was varied by changing the 
thickness of the rubber and by using natural 
rubber, Neoprene, and combinations of the two. 

III. Tests 
1. Release Gear 

The release gear shown in figure 7 was developed 
to drop the model in the nearly " strain-free con- 
dition" that should hold during free fall. This 
prevents the setting up of vibrations excited by 
the sudden removal of the dead-weight forces 
upon release of the model. It has been found 
that such vibrations may interfere seriously with 
the interpretation of the strains and accelerations 
recorded during drop tests of full-size airplanes. 

The release gear supports the model at several 
stations along the wing with forces that are ad- 
justed to be nearly in balance with the local dead- 
weight forces. Upon release, the supports are 
removed at an acceleration greater than gravity, 
thus leaving the model free to fall in its nearly 
strain-free condition. 

The support is applied to the model by eight 
pointed screws, B, at the ends of six arms, A, 
figure 7. The screws, B, were carefully adjusted 
until strain gages attached to the model near the 
root and near the first engine indicated the model 
was in a " strain-free condition". Strain readings 




Figure 6. Center section of idealized model. 

Journal of Research 




Figure 7. Release gear. 

corresponding to the "strain-free" condition were 
determined as the average of readings for the -\-g 
and the —g condition. The -\-g condition was 
obtained by holding the model above the center 
of gravity with the landing foot vertical below the 
center of gravity; the —g condition was obtained 
by turning the model through 180° about the 
wing axis to place the landing foot vertically 
above the center of gravity. 

The model was released by cutting an 0.04-in. 
steel wire attached to lever C, which held the 
release gear in the supporting position. The 
supports were swung away from under the model 
at an acceleration greater than g by heavy rubber 
bands, D, which rotated the main support rod, 
E, about ball bearings at the ends. A friction 
catch at F on rod E prevented the arms, A, from 
swinging back into the model. 

2. Instrumentation 

The instrumentation of the model included 
pick-ups for measuring bending moment near the 
fuselage and near the first outboard engine, accel- 
eration at the fuselage, force transmitted through 
the springs, and force transmitted through the 
damper. 

All quantities except acceleration were measured 
with pick-ups by using strain-sensitive wire. The 
signal of these pick-ups was amplified with the 



four-channel equipment shown at G in figure 7. 
It was then recorded in the six-channel recorder, 
H. The amplifying equipment included in each 
channel a 1,000 c/s carrier-current Wheatstone 
bridge and a band-pass filter to give nearly flat 
response up to 100 c/s; at 200 c/s the response 
was down about 20 percent. The recorder was 
equipped with galvanometers having a natural fre- 
quency of 430 c/s and a response that was flat up 
to 100 c/s.; the response rose about 20 percent 
at 200 c/s; it dropped off rapidly above 300 c/s. 

Acceleration was measured with a vacuum- 
tube acceleration pick-up developed at the Na- 
tional Bureau of Standards. The pick-up con- 
tains two separate flexibly mounted plates on 
opposite sides of a fixed cathode. The plates are 
deflected elastically by the accelerations, thereby 
increasing the current between one plate and the 
cathode and decreasing it between the other plate 
and the cathode. The total change in current is 
recorded through a Wheatstone bridge. The 
pick-up has a fundamental frequency of about 
800 c/s. The output was fed through a low-pass 
filter into the recorder, H, giving a flat response 
up to 200 c/s. 

The bending moments at various stations on the 
wings were measured by attaching pairs of wire- 
strain gages to the top and bottom of the wing 
and connecting them in opposition in a Wheat- 
stone bridge circuit so that the output was pro- 
portional to the extreme fiber-bending strain. 
In addition, gages in corresponding positions on 
the right and left halves of the wing were con- 
nected in series to average their output electri- 
cally. The over-all circuit from gages through 
Wheatstone bridge, amplifier, and oscillograph 
was calibrated statically before each set of drop 
tests by applying known bending moments to the 
wing and recording the output. 

The acceleration at the fuselage was measured 
with the vacuum-tube acceleration pick-up, A, 
figure 6. The pick-up and circuit were calibrated 
with an acceleration of 2 g applied before each 
set of drop tests by reversing the accelerometer 
in the earth's gravitational field and recording the 
output. 

The force transmitted through the spring was 
measured by attaching wire-strain gages, B, figure 
6, to each of the eight leaves of the spring and 
connecting them into a Wheatstone bridge circuit 
in such a way that their output was proportional 
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Figure 8. Drop tests from a height of about 0.7 in. 

Pressure relief ports open U turn in tests a to c. In test f, ports 1 to 5 open one-fourth turn, and ports 6 to 8 closed. A, fuselage acceleration; B, spring; 
C, damper; D, outboard bending; E, root bending. 
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to the total force transmitted by the spring. A 
static calibration was obtained before each set of 
tests by applying known forces and recording the 
output. 

The force transmitted through the damper was 
measured by a small pressure gage at the top of 
the cylinder chamber, figure 5. This pressure 
gage consisted of an aluminum alloy tube %-in. 
long, % in. in diameter, and of 0.0035-in. wall 
thickness. The tube had one end closed and the 
other end open to the fluid in the cylinder. On 
the outside of the tube 0.001-in. constantan wire 
was wound and cemented into place taking care 
to insulate the wire from the tube and to bring 
out firmly anchored lead wires. A static calibra- 
tion of the gage showed it to have a linear change 
in resistance with pressure up to pressures of 300 



lb/in. 2 , corresponding to forces of 500 lb trans- 
mitted through the damper. 

IV. Results 

The records showing the results of the drop 
tests from a height of about 0.7 in. are shown in 
figure 8. Six drops were made in which the land- 
ing conditions were varied by changing the softness 
of the rubber on which the model was dropped, 
the initial spring loading, and the opening of the 
pressure relief ports in the damper. The initial 
spring loading of 43 lb, for impact, (figure 8, /, 
corresponds to supporting the whole weight of 
the model on the springs. 

The shortest impact, lasting about 0.050 sec, is 
shown in figure 8, f, and the longest impact, last- 
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ing about 0.086 sec, is shown in figure 8, a. The 
ratio of impact time to the period of the model 
in its fundamental bending mode approximately 
covers the range of this ratio for large airplanes. 

The maximum observed bending moments in 
the wing at stations 1.5 in. and 14.5 in. from the 
root are given in table 1. In addition, table 1 
gives the maximum observed accelerations at the 
root, the duration of the impact, and the maximum 
impact force. The maximum impact force was 
obtained as the maximum value of the sum of 
the forces transmitted through the springs and 
the damper. 

Table 1. Maximum values of the applied impact force 
(sum of spring and damper forces), bending moments at 
1.5 -in. and ll+.5-in. stations, and accelerations at root 
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V. Analysis 

1. Normal Modes of Vibration 

The normal modes of vibration of the model 
wing in bending were computed by considering the 
mass to be concentrated at the root and at nine 
stations along each half-wing. The distributed 



mass of the portion of the wing between adjacent 
stations was distributed to those stations in in- 
verse proportion to the distance from each station 
to the center of gravity of the portion of wing 
being considered. The values are given in table 
2. Influence coefficients were then computed be- 
tween the mass points by treating each half-wing 
as a simple beam clamped at the root. The first 
three flexural modes of the wing in free-free vibra- 
tion were computed from these influence coeffi- 
cients and from the given mass distribution using 
a dynamic matrix and iteration procedure as ex- 
plained by Duncan and Collar [4]; the deflection 
at the root of the wing was obtained from the 
condition that the center of gravity of the wing 
must remain at rest for free-free vibration. In the 
case of the second and third mode, the iteration 
procedure had to be modified to prevent conver- 
gence to the fundamental mode. This was ac- 
complished by removing any proportion of lower 
modes present with the help of the orthogonality 
relations as outlined for propeller blades in [5]. 

The deflections in each mode r were normalized 
by dividing the deflection, y^, at a given sta- 
tion, i, by that at the tip station. The values are 
given in table 2. 

The normalized deflection, tf\ , were substi- 
tuted in the following formula to obtain the gen- 
eralized mass, M r , of the equivalent linear oscilla- 
tor in mode r (see [1]): 



i=0 



(i) 



where ra f is the mass at station i. The values are 
given in table 3, together with the frequencies and 
periods in each mode. 



Table 2. Distribution of mass and flexibility along half-iving ande shap of natural modes 
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Table 3. 



Natural period, generalized mass, and bending 
moment for each natural mode 
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Table 3 shows also the bending moments per 
unit- tip deflection in each mode at the two sta- 
tions, 1.5 and 14.5 in. from the root, at which 
measurements were made. The bending mo- 
ments were computed by considering the wing to 
be loaded at stations i=l, 2, . . ., 9 by transverse 
forces ufrrii rn (r) where co r is the frequency in the 
rth mode. 

2. Dynamic Response According to Theory of 
Biot and Bisplinghoff 

The dynamic response of the model wing was 
computed by using Biot and Bisplinghoff 's theory 
[1] for each of the 6 impacts shown in figure 9. 
These were obtained from the records of figure 8 
by adding the spring and damper forces shown as 
B and C and dividing by 2. 

(a) Bending Moment 

The maximum bending moments M M . 5 and 
M bM . 5 at stations 1.5 in. and 14.5 in. from the root 
were computed by adding up the maximum bend- 
ing moments (M bA , 5 ) r and (M 6t i 4 . 5 ) r for the first 3 
modes r=l, 2, 3. Two sets of values were com- 
puted, the first an upper limit corresponding to 
the envelope of dynamic response factors y r 
given in figure 13 of [1], and the second a closer 
approximation corresponding to the response factor 
for that impact in figure 12 of [1J that came 
closest to the actual impacts shown in figure 9. 

The bending moments in each mode r were 
obtained by multiplying the bending moments 
per unit tip deflection in table 3 by the tip deflec- 
tion y T Q r lM r <a r 2 f where Q r is the generalized force 
in mode r. The generalized force was computed 
as the product of the maximum observed impact 
force and the normalized deflection at the root 
in mode r. 




Time, seconds 



Figure 9. Total impact force per half wing for tests of 
figure 8. 

Table 4. Bending moments according to theory of Biot and 
Bisplinghoff, lb-in. 



Bending moments 



Drop 



(Mb, 1.5) 1 envelope . 
(Mb, 1.5) 2 envelope . 
(Mb, 1.5) 3 envelope - 



(Mb, 1.5) 1 best fit, 
(Mb, 1.5) 2 best fit. 
(Mb, 1.5) 3 best fit.. 



(Mb, 14.5) 1 envelope 

(Mb, 14.5) 2 envelope 

(Mb, 14.5) 3 envelope 

(Mb, 14.5) 1 best fit 

(A/6, 14.5) 2 best fit 

(Mb, 14.5) 3 best fit 

i 

(Mb, 1.5) envelope 

(Mb, 14.5) envelope 



(Mb, 1.5) best fit--. 
(Mb, 14.5) best fit.. 



-1.170 
-231 

-38 

-735 
-182 
-38 

-813 

-66 

3 

-511 

-52 

3 

-1,439 

-876 

-955 
-560 



c d 



-1,321 

-270 
-42 

-1,173 
-203 
-42 

-918 

-77 
3 

-815 

-58 

3 

-1,633 
-992 



-1,278 
-275 
-46 

-673 
-183 
-40 



-78 
4 

-467 

-52 

3 



— 1, £ 



■1,418 
-870 -516 



-1,182 
-256 
-43 

-650 
-175 
-37 

-821 

-73 

4 

-451 

-50 

3 

-1,481 



-862 



-1,550 
-349 
-63 



-242 

-49 

-1,077 

-99 

5 

-673 

-69 

4 

-1,962 
-1,171 

■1,260 
-738 



-1,460 
-331 



-226 
-47 

-1,014 

-94 

5 

-644 

-64 

4 

-1,851 
-1, 103 

-1. 201 
-704 



516 



Journal of Research 



The results of the computations are given in 
table 4. Examination of this table shows that 
the third flexural mode contributes less than 4 
percent to the bending moment at the two stations. 
In view of this rapid convergence, it would have 
been sufficient in this particular case to confine 
the analysis to the first two flexural modes. 

Comparison of the maximum bending moments 
computed, table 4, with those observed experi- 
mentally, table 1, shows that the computed values 
using the best fit impact curve are from 15 percent 
less to 51 percent more than those observed, and 
using the envelope impact curve are from 43 to 
137 percent more than those observed. 

(b) Acceleration at Root of Wing 

A dynamic response factor a for acceleration, 
analogous to the dynamic response factor y for 
deflections as given in figure 13 of [1], was deter- 
mined as follows: 

From equations (1-8) and (1-10) of [1], 



y^ T) ^ = Q M r -^ [ Qr(r) sin c r (t-r)dr, (2 



where 



#tip (r) (0 =acce l era ti° n in mode r at tip 

Q T (t)= generalized force in mode r as a 
function of time 
r = variable of integration 
t = tlme. 

If we denote by Q T the maximum value of Q r (t) and 
let p r (t) be a unit impact force defined by 



?#)=-£(}#), 



(3) 



we obtain the dynamic response factor a ir) (t) for 
acceleration as 



w !/ '' 



(0= 



2/tip 



in 



Qr/M r 



(t)=Pr(t)—(ti 



r J Pt(t) 



(4) 



The dynamic response factor a {r) (t) can be con- 
sidered as the ratio of the actual acceleration 
2/ti P (r) (0 a t the tip to the steady acceleration result- 
ing from applying the peak generalized force Q T 
to the generalized mass M T . 

The value of a {T) (t) was evaluated for impacts 
of triangular shape and sinusoidal shape by using 
the analyses presented in eq 16 to 19 of [6]. The 
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peak values a of a iT) (t) are plotted in figure 10 
as a function of the ratio of impact pulse period 
T to mode period T t =2-k\u t . The dynamic re- 
sponse factor a for acceleration decreases rapidly 
ns the ratio T/T r is increased above 0.8. For 
T/T T greater than 7, a is less than 0.20. As T/T T 
approaches zero, a approaches 1. 

The maximum acceleration y al the root of the 
wing was computed by adding up the maximum 
accelerations y Q (r) for the first four modes r=0, 1, 
2, 3, where r=0 corresponds to motion as a rigid 
body. Two sets of values were computed, as in 
the case of the bending moments, the first corre- 
sponding to the envelope of dynamic response 
factors a in figure 10 and the second corresponding 
to that one of the two shapes of the curve of 
impact force versus time, assumed for figure 10, 
that comes closest to the actual impact shown in 
figure 9. 

The accelerations y {r) in each mode were ob- 
tained from 



#o 



<r) = 



aQrVo 



(r) 



M r 



(5) 



where rj {) {r) is the normalized deflect ion at the root. 
For the rigid body mode r=0. 

y m = QJM , 



(5a) 



where (J {] is one-half of the maximum impact force 
applied, and M is one-half of the mass of the wing. 
The results of the computation are given in 
table 5. Examination of table 5 shows that the 
second and third flexural modes contribute less 
than 4 percent to the acceleration at the root. 
It would have been sufficient in this case to con- 

Table 5. Acceleration at root of wing according to theory 
of Biot and Bisplinghoff in. I sec 2 



Acceleration 


Drop 


a 


b 


c 


d 


e 


f 


jfo<°> 


1,720 

214 

52 

9 

180 
29 
3 

1,995 

1,932 


1,900 

257 

63 

12 

132 

43 

6 

2,232 

2,031 


1,810 

279 

74 

14 

78 
48 
4 

2,177 

1,940 


1,665 

259 

69 

14 

80 
62 
4 

2, 007 

1,813 


2, 172 

367 

105 

21 

201 
56 
11 

2,665 

2,440 


2,046 

357 
102 
20 

217 
41 
6 

2, 525 

2,310 


jro (^(envelope) .. 
yo (2) (envelope) . _ 
#o (3) (envelope) . _ 

fo a) (best fit). — 
jfo < 2 > (best fit).. .. 
jM 3 > (best fit).... 

jfo (envelope) 

£o (best fit) 
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Figure 10. Dynamic response factor for acceleration. 
Triangular pulse; , sinusoidal pulse; , envelope. 



fine the analysis to the rigid body motion and the 
fundamental flexural mode. 

Comparison of the maximum computed accel- 
erations at the root, table 5, with those observed 
experimentally, table 1, shows that the computed 
values using the envelope curve are from 16 to 
56 percent more than those observed and, using 
the best fit curve, are from 8 to 43 percent more 
than those observed. 

3. More Exact Analysis, Including Effect of Phase 
Differences 

A basic assumption in the theory of [1] is that 
the maximum response of the wing in each mode 
may be added without regard for phase. The 
resulting error was computed by making a more 
exact analysis for the impact-time relation shown 
in figure 9, d. 

The tip deflection in mode r(r=l, 2, 3) was 
computed from eq (I — 10) of [1] 



y^w^r Jo^ r(r) sin w ' ( *~ r)rfr > 



(6) 



where Q r (t) is the generalized force in mode r 
computed by multiplying the impact force-time 



relation figure 9, d by the normalized deflection at 
the root, r$ . The integrations were carried 
out numerically. The resulting tip deflections in 
each of the first three flexural modes are shown in 
figure 11, a. 

The bending moments M b , U5 at 1.5 in. from 
the root were computed by multiplying the tip 
deflections in each mode by the corresponding 
factor in table 3. The resulting bending moments 
in each of the first three modes and their sum, as 
well as the observed bending moment, are shown 
in figure 11, b. Figure 11, c shows corresponding 
values at 14.5 in. from the root. 

The acceleration at the root in the rigid body 
mode r=0 was computed from eq 5a after replac- 
ing Q by Qo(t). The acceleration in the flexural 
modes r=l, 2, 3 was obtained from 



2/o 



(')- 



"2/tip Vo 



[^-"rW"]*"' (7) 



where y^l is the tip deflection given in figure 
11, a. The results of the computation are given 
in figure 1 1 , d, together with the observed accelera- 
tion at the root. 

Examination of figures 11, b to 11, d shows 
that the observed maximum bending moments 
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Figure 1 1. Detailed analysis for the impact force-time re- 
lation in figure 9,d. 

and accelerations differed less than 20 percent 
from the computed values, when the phase differ- 
ences were taken into account and when the 
actual impact force-time curve was used. 

Comparison of the response in the individual 
modes in figure 1 1 shows that the error caused by 
neglecting phase differences was negligible (less 
than 1%) in the case of bending moments. In 
the case of acceleration at the root, the error was 
about 23 percent. 

The good agreement observed for bending 
moment is primarily the result of the following 
two factors. 

In aircraft structures (or models simulating air- 
craft), the periods of natural vibration T r are as 
short or shorter than the duration of the impact 
T. In such a case, a maximum bending moment 
and deflection response in all the modes occurs in 
phase at a time a little later than the maximum 
applied impact force. The response in the differ- 
ent modes drops off rapidly as the number of the 

Landing Impact 



mode increases. Hence the contribution of the 
higher modes becomes negligible. 

The phase difference can be important in the 
case of acceleration, as the acceleration in the 
first flexural mode may be, and in this example 
was, opposed in phase to that of the rigid body 
motion at the peak. Neglect of phase in the case 
of acceleration near the tip of the wing would 
probably lead to results so inaccurate as to be 
useless. 

VI. Summary 

An experimental verification of Biot and Bis- 
plinghoff's analysis of landing impact is presented. 
An airplane model was built having a distribution 
of mass and of flexural rigidity along the wing ap- 
proximately proportional to that for a four-engine 
military airplane. The four engine masses were 
mounted symmetrically on the wing so as to 
excite flexural vibrations without torsion when the 
model was dropped vertically to receive a landing 
impact below the center of gravity. The model 
contained an alighting gear with means for adjust- 
ing the time history of the impact force acting on 
the wings. Measurements were made of impact 
force, bending moments at two stations, and root 
acceleration for six landing conditions. 

The observed maximum bending moments were 
compared with those computed by Biot and Bis- 
plinghoff's method. The computed bending mo- 
ment, using a response factor given by an envelope 
curve due to Biot and Bisplinghoff, was 43 to 137 
percent larger than the observed bending moment. 
Use of an approximation to the actual impact-time 
relation gave bending moments from 15 percent 
less to 51 percent greater than the observed bend- 
ing moments. 

A more exact analysis, taking account of phase 
differences was made for one of the actual impact 
force-time curves. This gave maximum bending 
moments that differed less than 20 percent from 
the observed bending moments. Comparison of 
the computed response in the individual modes 
showed that the neglect of phase differences be- 
tween the different modes, in accordance with 
Biot and Bisplinghoff ' s analysis, would have 
caused an error in maximum bending moment of 
less than 1 percent. 

The observed maximum accelerations at the 
wing root were compared with those computed 
from an extension of Biot and Bisplinghoff's anal- 
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ysis, which involved the computation of response 
factor curves for acceleration. The computed 
acceleration, using a response factor given by the 
envelope to these curves, was 16 to 56 percent 
greater than the measured acceleration. Use of 
the most closely fitting impact curve reduced this 
error to 8 to 43 percent. It should be noted in 
this connection that larger differences may be ex- 
pected at the wing tips, since these are more 
affected by flexural vibrations. 

The more exact analysis for one of the impacts 
showed that the neglect of phase differences caused 
an error in computed maximum accelerations at 
the root of about 23 percent. This error would 
probably be much greater for acceleration at the 
wing tip. 

The authors thank their colleagues Samuel Levy, 
for the analysis reported in this paper, and Eileen 



D. Segal, for the numerical computations of dy- 
namic response. 
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